On the density matrix of nonequilibrium steady-state statistical mechanics 
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This paper derives a density matrix of the steady-state statistical mechanics compatible with the 
steady-state thermodynamics proposed by Oono and Paniconi [Prog. Theor. Phys. Suppl. 130, 29 
(1998)]. To this end, we adopt three plausible basic assumptions for uniform steady states: (i) 
equivalence between any two subsystems of the total, (ii) statistical independence between any two 
subsystems, and (iii) additivity of energy. With a suitable definition of energy, it is then shown that 
uniform steady states driven by mechanical forces may be described by the Gibbs distribution. 
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Constructing thermodynamics and statistical mechan- 
ics far from equilibrium is undoubtedly a major goal yet 
to be achieved in modern condensed matter physics. One 
strategy tovi^ard this may be to seek a way to extend 
the well-established equilibrium framework to nonequi- 
librium systems. Numerous efforts have been made along 
this line|, I I, I, I, |, 0, I . However, most of 
them starts from the local-equilibrium hypothesis which 
may not be justified for systems far from equilibrium. 
Recently, Oono and Paniconi Q presented a new ap- 
proach restricting their attention to nonequilibrium time- 
independent states. A key ingredient lies in the removal 
of "house keeping heat rate" Qhk which is generated in 
the system as dissipation to be carried away eventually 
by some microscopic degrees of freedom. Using the "ex- 
cess heat rate" Qex defined by subtracting Qhk and con- 
necting points of steady "state space" in a well-defined 
way, they have constructed a thermodynamic framework 
named "steady-state thermodynamics" (SST) which is 
quite analogous to the equilibrium one. Indeed, the en- 
ergy and the Helmholtz free energy are defined there in 
the same way as equilibrium thermodynamics by merely 
introducing extra extensive variables characteristic of the 
relevant steady state. 

Naturally, an amount of their effort was directed to- 
wards an unambiguous definition of entropy which serves 
as the potential for nonequilibrium systems to determine 
their evolution and stability criteria. It should be noted 
that they are completely different from the variational 
principle on entropy production rate established by On- 
sager near equilibrium [p^ and taken up by Glansdorff 
and Prigoginc|^ as a candidate for the general criteria 
in nonequilibrium systems. Indeed, entropy and entropy 
production rate are different in dimension so that the 
two variational principles cannot be compatible with each 
other. It is intuitively plausible to expect that, once the 
steady dissipation is thrown away into "hidden degrees 
of freedom" |pT| , we may have well-defined energy and en- 
tropy which play the same role as the equilibrium thermal 
physics. 

Efforts have been made thereafter to test SST based on 
microscopic stochastic models |p2|, p^ , p^ . However, 
those stochastic models are so designed as to settle down 
in equilibrium, it is not clear whether the results from 



them are really relevant to far-from-equilibrium states, 
as recently pointed out by Sasa and Tasaki[p^. 

We here take an alternative approach to seek for a 
statistical density matrix which is compatible with SST 
of Oono and Paniconi ||^. There are a couple of mer- 
its in this approach. First, once such a density matrix 
is obtained, we may perform microscopic calculations on 
nonequilibrium steady states which could be tested by ex- 
periments. Second, this statistical mechanics may make 
clearer the concepts of SST such as "state space," "en- 
ergy," and "entropy." Hopefully the attempts will also be 
helpful to clarify the conditions under which SST holds. 

We restrict ourselves to the cases of mechanical ex- 
ternal perturbations where an unambiguous treatment is 
possible. As is well known, the thermal perturbations 
are difficult to handle even within the linear-response 
regime[0, so that we leave those cases for a future con- 
sideration. 

We specifically consider electrons in a metal with im- 
purities which is embedded in a uniform electric field E; 
this electric field sustain uniform current J in the sys- 
tem. First, let us review a conventional way to describe 
the system: It is convenient to express the electric field 
in terms of the time-dependent vector potential A.{t) as 
E = — ^"^1*'' ■ The corresponding Hamiltonian is given in 
units of ^ = 1 by 



n{t) = U\v) 



2m 



-Kinp(r)-/i >V(r) dr 



(1) 



with e (< 0) the electron charge, m the electron mass, 
Vimp the impurity potential, and /i the chemical po- 
tential. The spin degrees of freedom and the electron- 
electron interaction are suppressed for simplicity. A gen- 
eral advantage of using time-dependent Hamiltonians for 
nonequilibrium states is that we can definitely identify 
the input power as (^^p^)- Here, electrons are acceler- 
ated along E to acquire the input power, but then scat- 
tered by Vimp resulting in their momentum relaxation. 
Finally, the energy relaxation occurs, i.e., the electrons 
interact with some other microscopic degrees of freedom 
such as phonons, and the extra energy accumulated in 
the electrons are eventually carried away from the sys- 
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tern. However, this last stage cannot be described by 
the density matrix p{t) obtained by solving the quantum 
Liouville equation 0, p7[|: 



Indeed, its formal solution is given by 
pit) = Sit)poS\t), 



(2) 



(3) 



where po is the equilibrium density matrix at t = Iq when 
E is switched on, and S{t) is defined by 



S{t) = Texp 



i / n{t')dt' 



to 



(4) 



Equation (^) forms a basic starting point for the 
linear-response theory and the fluctuation-dissipation 
theorem^, Even beyond the linear-response regime, 
it provides a general framework to perform nonlinear 
nonequilibrium calculations However, since p(t) 
develops purely mechanically from equilibrium po with 
no thermal contact, energy relaxation processes are ab- 
sent in Eq. (^). Thus, exact calculations based on Eq. 

predict that the energy of the electron system grows 
towards infinity as time goes by. It might be possible to 
add terms responsible for the energy relaxation to Eq. (j^) 
Q , but there seems to be an ambiguity in this procedure 
especially in the nonlinear region. Thus, the quantum 
Liouville equation may not be a good starting point to 
obtain the correct steady-state density matrix far from 
equilibrium. Even more difficult will be to perform calcu- 
lations of including microscopic degrees of freedom which 
transport the extra electron energy out of the system, be- 
cause whether the combined system settles down into a 
steady state by itself is not entirely clear. 

To seek for the density matrix, we here proceed differ- 
ently based on an argument which assumes the existence 
of a uniform steady state from the beginning and which 
fully relies upon the uniformity of the system. In the 
end, the equilibrium density matrix can be derived from: 

(i) equivalence between any two subsystems of the total, 

(ii) statistical independence between any two subsystems, 
and (iii) additivity of energy. These are essential ingre- 
dients of the concept of entropy. Since (i) and (ii) are 
also characteristic of any uniform steady states, we may 
expect to have a well-defined entropy by a suitable def- 
inition of energy for each subsystem. An extension to 
nonuniform steady states may be performed through the 
usual procedure of dividing the total system into small 
uniform cells ||. 

Consider a subsystem a of the total where steady cur- 
rent J is present (see Fig. 1). In this subsystem, there 
is equal amount of energy input and dissipation per unit 
time, the latter being thrown out of the subsystem by 
some microscopic degrees of freedom not considered ex- 
plicitly. Following the philosophy of SST, we disregard 




FIG. 1: The total system with a uniform current and its sub- 
systems. 



this energy flow, include the carrier of the dissipations in 
"hidden degrees of freedom" , and consider the fluc- 
tuations caused by the energy flow and those through 
the wall of the subsystems together. Now, the system 
is characterized by the uniform static current J, whose 
"energy" may be calculated by the static "Hamiltonian:" 



(-zV-po) 
2m 



+ Vinip(r) - p 



ip{r) dr , 



(5) 



where po is a variable conjugate to J. From now on 
we can follow the argument of equilibrium statistical me- 
chanics to derive the density matrix p^"-^ of the subsystem 
a: Since it is time-independent by assumption, p^'^^ com- 
mutes with Ti.'-°'\ so that it is a function of the eigen- 
value si)^^ of H*^"). From (ii) the statistical indepen- 
dence, = p'^°-^ p^^\ and (iii) additivity of energy, 

g{a+b) ^ g(a) _^ g{b)^ j^- follows thatpl 



(6) 



where a^"^^ and (3 are some constants. Normalizing p^"', 
we have 



E 



(7) 



where higher probability for lower £ means (3 > Q. Thus, 
we have reached a candidate for the steady state density 
matrix in the presence of a uniform current J. 

We now see that Eq. (0) corresponds to the state 
of maximum entropy, following exactly the argument 
of equilibrium case po|| : Entropy is defined by S*'"' = 
— Trpf") In = — (inp^"') and calculated by using Eqs. 
d) and d) as 

= + /3((^('')) - p'(iV('')) - Po • (J("))) , (8) 

where is the energy of the subsystem, (J'^"-') = 

-(aW^'^V^Po), and p' = p + We also notice that 
-(lnp('^)) = - In ((£)). On the other hand, Eq. (g) 
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tells us that the probability is the same for all degener- 
ate states and equal in the thermodynamic limit to the 
inverse of the number of states around Hence 
Eq. corresponds to the entropy maximum for the en- 
ergy the electron number (A^'^"^), and the current 
Finally, we may perform the Legendre transfor- 
mation G^") = (H^")) + po • (J^°^) to obtain the desired 
free energy G^"^ as a function of /?, fi', and (J^"^^). 

The above consideration can be applied to any uni- 
form steady states driven by mechanical forces. Indeed, 
we only have to introduce in the Hamiltonian a variable 
[such as Po in Eq. (^)] which is conjugate to the expec- 
tation value (such as J in the above consideration) char- 
acterizing the steady state. 

Several comments are in order. First, the considera- 
tion here puts aside completely the driving force E and 
the corresponding dissipations. Thus, this density ma- 
trix cannot say anything about the E-J relation, but only 
identifies the steady state for a given J. However, this 
E-J relation may be obtained by connecting the results 
through the linear-response calculations for each J. Sec- 
ond, the above argument starts from three plausible ba- 
sic assumptions on uniform steady states: (i) equivalence 
between any two subsystems of the total, (ii) statistical 
independence between any two subsystems, and (iii) ad- 
ditivity of energy. Hence we may expect that the derived 
Gibbs distribution (|^) is the correct distribution for those 
steady states. Indeed, there are several stochastic mod- 
els which realize the Gibbs distribution in nonequilibrium 
steady states. However, it should be noted that 

the quantity /? here includes the fluctuation of energy due 
to dissipations as well as that through the walls among 
subsystems. If the former contribution is negligible, we 
may expect that (3 is the same as the equilibrium value. 
Experimentally, there remains a basic problem on how 
to measure the value of (3 for the steady states. Third, 
the above consideration have surely made the concept 
of "state space" in SST clearer. The postulated second- 
law of SST may be proved now by using Eq. (|^) and 
the Jarzynski identity[Q. Fourth, the subsystem con- 
sidered here may contain structures small compared with 
the subsystem size such as the vortex lattice structures. 

There are many interesting nonequilibrium phenomena 
found in semiconductors. [25| Also, StoU et al. found steps 
in the I-V characteristics in the vortex state of supercon- 
ducting Nd2-2:Ce2,CuOj^; this first-order-like transition 
may be described by the nonequilibrium free energy de- 
rived above, as increasing E causes the density-of-states 



(DOS) change in the electron system. 

I am grateful for Koji Nemoto for an enlightening dis- 
cussion. 
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